Stable nondegenerate optical parametric oscillation at degenerate 

frequencies in Na:KTP 



Sheng Feng and Olivier Pfister* 



m 
o 
o 

(N 

Ph. 
< 



(N 
> 



(N 
O 

:^ 



Department of Physics, University of Virginia, 
382 McCormick Road, Charlottesville, VA 22904-4714, USA 

February 1, 2008 



Abstract 

We report the realization of a light source specifi- 
cally designed for the generation of bright continuous- 
variable entangled beams and for Heisenberg-limited 
inteferometry. The source is a nondegenerate, single- 
mode, continuous-wave optical parametric oscillator 
in Na:KTP, operated at frequency degeneracy and 
just above threshold, which is also of interest for the 
study of critical fluctuations at the transition point. 
The residual frequency-difference jitter is ± 150 kHz 
for a 3 MHz cold cavity half-width at half maximum. 
We observe 4 dB of photon-number-difference squeez- 
ing at 200 kHz. The Na:KTP crystal is noncritically 
phase-matched for a 532 nm pump and polarization 
crosstalk is therefore practically nonexistent. 

1 Introduction 

Optical parametric oscillators (OPO's) are well 
known sources of nonclassical light. Type-I and type- 
H OPO's have been used to generate single- mode 
and two-mode [3 EI squeezed states and to imple- 
ment the Einstein-Podolsky- Rosen (EPR) paradox |1| 
experimentally as proposed in [Sj The lat- 
ter has become increasingly important as continuous 
variables (CV) have been proposed and used to im- 
plement quantum information protocols [7||S| such as 
quantum teleportation 9 and quantum dense coding 
[TU| . This paper is focused on CV and OPO's, but 
one should not forget, of course, that single-photon 
coincidence detection coupled with spontaneous para- 
metric downconversion ^Jj has played a major role 
in the violation of Bell's inequalities (practically 
unfeasible with CV), and the realization, for exam- 
ple, of teleportation ^3], multipartite entanglement 
|14j , and very recently quantum computer gates |15j . 
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Going back to CV, the generation of CV-EPR 
states requires simultaneous squeezing of two com- 
patible joint operators made from noncompatible lo- 
cal operators ^ [7], or physical variables, such as 
number sum (difference) and phase difference (sum) 
(iV±,0ip), or amplitude-quadrature sum (difference) 
and phase-quadrature difference (sum) (X±,Pzp). 
Reid and Drummond proposed to use above- and 
below-threshold OPO's in order to squeeze N- and 
4>+ IHl In the first experimental implementa- 
tion of the EPR paradox a nondegenerate OPO 
below threshold was used. An interesting variant 
also demonstrated was to use two degenerate OPO's 
below threshold, each squeezing a complementary 
quadrature, and to overlap their fields using a bal- 
anced beam splitter jH] . As pointed out by Lugiato et 
al., the unitary transformation describing the beam 
splitter transforms the nondegenerate and degener- 
ate Hamiltonians into one another |16| . Several re- 
lated entanglement experiments have been recently 
reported |17[ I18L and one ought also to mention 
the use of third-order nonlinearity in fibers to gen- 
erate bright EPR beams without the use of OPO's 

It is interesting to note that, despite the fact that 
the initial proposal for experimentally realizing the 
EPR paradox was to use an OPO above threshold 
[5], this particular implementation has not yet been 
carried out. Rather, in all EPR implementations to 
date El El [17, 18, 19., OPO's have been used ex- 
clusively below threshold, as parametric amplifiers of 
either vacuum or bright modes. The reason for this 
choice is that the proposal of Ref. jS] is very challeng- 
ing. It requires either a triply resonant OPO with 
higher finesse for the pump than for the signal and 
idler, which presents technical difficulties, or a doubly 
resonant resonant OPO just above threshold. More- 
over, in both cases, the OPO must be emitting close- 



1 



to-degenerate signal and idler frequencies, which is 
difficult to achieve because the well-known clustering 
effect in a doubly resonant type-II OPO puts strin- 
gent stabilization requirements on the crystal tem- 
perature, the cavity length, and the pump frequency. 
One has thus to address two experimental challenges, 
using a doubly resonant OPO: to achieve stable os- 
cillation on the frequency degenerate mode and do so 
close to the threshold, where the output power is the 
least stable. 

We have solved these two issues and present, in 
this paper, the experimental realization of a stable 
type-II OPO, emitting at frequency degeneracy while 
pumped only a few percent above threshold. Note 
that stable twin beams have been achieved before, in 
particular in the record-breaking number-difference 
squeezing experiments by the group of Fabre and Gi- 
acobino and in the more recent work of several 
other groups [21], but our present work is, to the 
best of our knowledge, the first time that stable twin 
beams are obtained at precise frequency degeneracy 
and just above threshold. 

Besides bright EPR beam generation, such a source 
is also extremely interesting for the study of critical 
fluctuations at the transition point, including possi- 
ble generation of macroscopic entanglement and, 
last but not least, for Heisenberg-limited interferom- 
etry (HLI). Recall that, if N particles are sent into 
an interferometer and subsequently detected or mea- 
sured, the usual phase sensitivity limit is the input 
beam splitter's shot noise limit N~^/^ [21], whereas 
the ultimate limit is the Heisenberg limit j^HI 
126) ■ This is a general statement for any boson field. 
Examples of beam splitters are a half-reflecting mir- 
ror in photon optics and a 7r/2 resonant laser pulse 
in atom optics. Reaching the Heisenberg limit re- 
quires suppressing vacuum fluctuations at the input 
beam splitter of the interferometer [25] . which was 
experimentally demonstrated using vacuum squeez- 
ing [2Z|. Holland and Burnett [23 a-lso showed that 
one can reach the Heisenberg limit by using a source 
of indistinguishable twin modes, i.e. an input density 
matrix of the form p = ^ pnm \ nn) ( m to | . This 
was demonstrated using two single-mode amplitude- 
squeezed beams in Ref. j20| . and also a pair of trapped 
ions [2Hi- Progress has also been made towards real- 
izing HLI with Bose-Einstein-condensate Fock states 
|29j . It is interesting to note, however, that the 
common-mode statistics {pnm)n.m do not play any 
role in the phase sensitivity |30| . hence generating 
twin Fock states e.g. \nn) is not necessary. There- 
fore, an OPO emitting intense frequency-degenerate 
twin beams is an ideal candidate system for HLI, as it 
would bring, besides the required iV_ squeezing, the 



added benefit of a a narrow linewidth CW source, 
suitable to perform ultra-precise measurements such 
as gravitational-wave detection. Note also that, since 
HLI requires large photon numbers, it presents one 
fewer challenge than bright EPR state generation. 

In the next section, we describe our experimental 
setup. We subsequently detail our experimental char- 
acterization of the source. 

2 Experimental setup 

2.1 OPO crystal and cavity 

The Na:KTP crystal noncritically phase-matches 
type-II degenerate parametric downconversion from 
532 to 1064 nm wavelengths along the X axis at 
room temperature. Such a configuration practically 
annihilates polarization crosstalk while allowing the 
convenient use of highly stable Nd:YAG lasers. We 
use a 3 X 3 X 10 mm"^ prototype Na:KTP crystal, 
fabricated by Crystal Associates. The crystal is cut 
along its X axis, and AR-coated at 1064 nm and 532 
nm. It is placed inside an oven made of Oxygen- 
free, high-conductivity Copper, whose temperature is 
controlled to less than a millidegree by way of a servo 
loop. We used two kinds of loop filters with equal suc- 
cess, one home-made and one commercial from Wave- 
length Electronics. Tight temperature control of the 
nonlinear crystal is paramount in a type-II OPO be- 
cause of the different temperature dependences of the 
indices of refraction for the two cross-polarized reso- 
nant signal and idler waves. 

The OPO has a standing-wave cavity formed by 
two plano-concave mirrors that have a radius of cur- 
vature of 5 cm and are separated by approximately 
10.2 cm, thereby realizing a concentric resonator. 
One of the mirrors is mounted on a piezoelectric 
transducer (PZT). The cavity waist size is in between 
the Boyd-Ashkin confocal- focusing value Wo{BA) = 
y^£X/2TTn = 31 pm and the Boyd-Kleinman opti- 
mum Wo{BK) ~ 0.59 ?i'o(-B^) = 18 pvn 31 , which is 
confirmed by OPO threshold measurements. The in- 
terferometer structure is built in super-Invar, which 
confers outstanding mechanical stability along with 
no visible thermal expansion effects. This makes for 
an ultrastable optical resonator that never needs to 
be realigned, the OPO threshold staying at its initial 
level for months on. 

The concentric geometry is more delicate to work 
with than a plano-concave cavity (as used for exam- 
ple in hemilitic OPO's) since its spherical symmetry 
allows for many different cavity axes and does not 
constrain the alignment very strongly. This is, how- 
ever, a conscious design choice as it allows us to eas- 



2 



ily find a cavity mode that coincides with the crys- 
tal's X axis so as to completely eliminate walkoff. 
The OPO alignment procedure consists of aligning 
the resonator without the crystal first, and then in- 
serting the crystal and orienting it between crossed 
polarizers in order to align the X axis with the cav- 
ity axis. The cavity is then used as a resonant fre- 
quency doubler in order to optimize the crystal tem- 
perature for maximum SHG efficiency, which corre- 
sponds to the lowest OPO threshold. The SHG pump 
laser at 1064 nm (Lightwave Electronics 126-700- 
1064) is, once frequency-doubled, heterodyned with 
the OPO pump laser (Lightwave Electronics 142) and 
the resulting beat note at their frequency difference 
is tuned to zero. The reflectivities of the OPO mir- 
rors are i?(1064 nm) = 0.999 and R{532 nm) = 0.6 
for the input coupler, and i?(1064 nm) = 0.990 and 
i?(532 nm) = 0.1 for the output coupler. The OPO 
is therefore essentially doubly resonant for the signal 
and the idler. 

2.2 OPO stabilization 

In order to ensure single-mode oscillation within the 
dense cluster structure (see below), we use stan- 
dard laser stabilization techniques j32| . which per- 
mit cavity-length control at the femtometer level ■ 
The OPO cavity is electronically stabilized by two 
servo loops, one for temperature and one for cav- 
ity length (Fig. The cavity-length error signal is 



532 rm 
pjmp laser 








Na:KT P^^ PBS 

T InGaAs PIN 



(D 



PIN 

<, / pholodiodes 



Loop filter 



Sp^njm 
analyzer 



Figure 1: Experimental setup. DM: Dichroic mirror 
(reflects 532 nm; transmits 1064 nm). 

picked off the small leak through the HR input cou- 
pler. A sinusoidal voltage at 25 kHz excites the first 
mechanical resonance of the PZT-mirror system and 
serves to modulate the OPO frequency. The weak 
leak signal is demodulated by a lock-in amplifier, pro- 
cessed by a loop filter, and finally applied to the OPO 
PZT. The good signal-to-noise ratio of this method 
allows us to operate with the OPO a few percent 



above threshold, which is of interest for quantum op- 
tics with bright beams as it will allow us to easily 
perform homodyne squeezing measurements with a 
stronger local oscillator, something believed difficult 
|17| because of the current power-withstanding limi- 
tations of high-efficiency photodiodes. 

2.3 Detection 

The OPO signal and idler beams are coUimated by 
an AR-coated lens and then traverse a variable beam 
splitter comprised of a half-wave plate followed by a 
low loss polarizing beam splitter {Tp, Rg > 0.99). We 
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Figure 2: (a): Intensity difference spectrum of twin 
beams through a balanced beam splitter (half-wave 
plate @ 7r/8 rad). The 4 MHz signal-idler beat note 
is clearly visible, along with its second harmonic 
(preamplifier distortion), (b): Intensity difference 
spectrum of twin beams (half- wave plate @ rad). 
The residual 4 MHz beat note is 52 dB below its 
maximum level of fig. [21 a) . (No intensity squeezing 
is visible here because the OPO beams are attenu- 
ated before detection.) Resolution bandwidth: 100 
kHz. Single sweeps. 

determine the reference angle of the wave plate axes 
a = by minimizing the beat note on the detectors. 
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which clearly separates the twin beams. For other 
angles relative to this one, the beams are thus mixed 
in the equivalent way to a beam splitter with reflec- 
tivity R = cos^(2a). When a = |, the setup is equiv- 
alent to a balanced beam splitter. This arrangement 
also allows us to check the presence of polarization 
crosstalk, i.e. the existence of a residual beat note for 
a — 0. As seen in Fig. El this effect is negligible. 

The detection is performed by two high effi- 
ciency (> 95%) InGaAs photodiodes (JDS Uniphase 
ETX500T) whose photocurrents are subtracted and 
sent into an ultralow-noise transimpedance preampli- 
fier. We measure the common-mode rejection ratio 
to be 42 dB @ 100 kHz, which enables us to have a 
shot-noise limited detection bandwidth as low as 200 
kHz. 



3 Experimental results 

3.1 Intensity-difference squeezing 

Figure 121 displays the measured photocurrent spec- 
trum for different cases, with a maximum measured 
squeezing of 4 dB @ 200 kHz. This detection band- 
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Figure 3: OPO intensity-difference squeezing spec- 
trum, (a): shot noise level — two-detector differ- 
ence signal in beam splitter mode (half-wave plate @ 
7r/8 rad). (b): number-difference squeezing — two- 
detector OPO difference signal (half-wave plate @ 
rad). (c): electronic detection noise (no light). Reso- 
lution bandwidth: 10 kHz. 100-sweep averages. 

width is quite low and unusual for squeezing ex- 
periments (but see also Ref. jSH). Reaching even 
lower quantum-limited detection bandwidths (1 kHz 
or less) would be of tremendous interest for HLI ap- 
plied to gravitational-wave detection. The squeezing 
spectrum has a maximum of 4 dB and an inverted 



Lorentzian shape of width the cold cavity linewidth 
j2|, which is 3 MHz here. Several straightforward 
measures can be taken to improve on this prelim- 
inary result, such as removing the photodetectors' 
uncoated windows, using a shorter crystal, and us- 
ing ultralow-loss OPO mirrors |3S1. The focus of this 
paper, however, is not as much on squeezing as on 
stability. 

3.2 Frequency control 

It is well-known that, in a doubly resonant type-II 
OPO, the double resonance condition for the cross- 
polarized signal and idler inside a common cavity 
containing a birefringcnt nonlinear crystal leads to 
a densely clustered mode structure This leads 

to a several hundred-fold increase of the constraints 
for length stability of the cavity, and demands that 
an electronic servo-lock loop be used. Theoretical 
analyses of the type-H OPO have been carried out in 
Refs. jSZl |3H1 so we just outline the relevant points 
here. The double resonance condition for the sig- 
nal and the idler is L -I- ris,i^ = Ps,iAs,i/2, L being 
the cavity length outside the crystal, £ the crystal 
length, Us^i the refraction indices of the signal and 
idler waves, Xg^i their respective wavelengths, and 
Ps^i their respective mode numbers (positive inte- 
gers). The equivalent resonance conditions for the 
optical frequencies Vg^i = c/Xs,i are iys,i — Ps,i^s,i = 
Ps,ic/[2{L -\- ris^ii)], where A^^^ are the free spectral 
ranges (FSR) of the signal and the idler. The reso- 
nance frequency change for a small air path change, 
which corresponds to actuating one of the mirrors 
with a PZT, is therefore 
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From this equation, one finds the well-known mini- 
mum length variation SL — ±A/2 for which a mode 
hop (5p = ^1, 5v — 0) occurs, in the case of a 
singly-resonant OPO or laser. In a doubly-resonant 
OPO, the additional constraint that both frequen- 
cies sum up to the pump frequency Vp =^ Vg + vi 
modifies the length tuning of the OPO cavity in a 
dramatic way, as the OPO will not lase for all com- 
binations of the signal and idler modes. The result is 
a tightly packed cluster structure. For a stable and 
narrow pump frequency, which is the case here, one 
has d(ys^v;l) I dL = 0. The length displacements cor- 
responding to mode hops are then given by, to first 
order in the birefringence Sn = Ug — rii, 
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where Xp is the pump wavelength, n — (jig + ni)/2, 
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Figure 4: OPO beat note spectrum, 5 MHz span. 
The narrow peaks are single spectrum analyzer 
sweeps (14 ms), and the broad features are the corre- 
sponding "max hold" traces for one minute or more, 
(a): Free-running OPO 4 times above threshold. The 
HWHM of the narrow peak is 50 kHz and the drift 
range is 2 MHz. Free-running oscillation stays single- 
mode without servo because the OPO mode HWHM 
is of the order of 10 MHz (the cluster modes overlap. 
See text), (b): Frequency- locked OPO four times 
above threshold. The narrow HWHM is 60 kHz and 
the drift range is 310 kHz, only limited by the signal- 
to-noise ration of the servo error signal. Note that 
the fast beat-note linewidth is limited by the 30 kHz 
resolution bandwidth of the analyzer. 

and Sps,i are integers. Mode clusters are labeled by 
{Sps + Spi) and separated by half a pump wavelength. 
Inside a given cluster, the modes are labeled by {Sps — 
Spi) and separated by {Sps — —Spi — ±1) 



5n 



. + L/t 2 



I « Ap, 



(3) 



the signal-idler frequency difference or beat note i^b — 
i>s ~ i^i differing, between two consecutive modes, by 



the sum of the FSR's. In the case of our experiment, 
n = 1.8, 6n — 0.09, L — 9.2 cm, and Xp — 532 nm, 
give clusters every 266 nm with mode hops 



Xp 
244 



2.2nm. 



(4) 



We now derive the tuning coefficients of the OPO 
sum and difference frequencies iy± = Vs ^ i^i with 
respect to the pump frequency Vp, cavity length L, 
crystal temperature T, and voltage across the crystal 
V (electro-optic tuning). Using Eq. ^ and analogs 
and assuming — and no mode hops, one finds 



( )t.v,., 
(dj^X 

\ ^^P J L,T,V 



(-5.11,-0.02) MHz/nm; (5) 
(-2.12,0.24) MHz/mK; (6) 
(1.34,0.59) MHz/V; (7) 
(1,0). (8) 



The next question is: what are the sources of noise? 
In the experiment, the pump frequency noise is < 10 
kHz/ms with a slow thermal drift < 10 MHz/min; the 
temperature servo gives ST < 1 mK in a 1 Hz band- 
width; and the cavity servo gives SL = 0.01 nm in a 
3 kHz bandwidth. No electro-optic servo is applied 
as the noise on V is assumed negligible. The results 
of these stabilization efforts are displayed on fig. 0] 
which shows recordings of the beat note at two very 
different time scales, with the temperature controlled 
and a free-running and locked cavity. The measured 
fast linewidths are of the order of 50 kHz for a 30 
kHz resolution bandwidth, and are thus instrument- 
limited. The locked residual jitter linewidth is 310 
kHz, down from the free-running range of 2 MHz, 
and the jitter frequency is 100-500 Hz. This 310 kHz 
residual linewidth on the beat note is larger than ex- 
pected from the length and temperature tuning co- 
efficients calculated above, and is too fast for being 
imputed to residual temperature fluctuations. We de- 
duce that the beat note stability is degraded by other 
factors. 

In fact, we have observed that the crystal is slightly 
inhomogeneous, likely due to a dopant concentration 
gradient. As a result, moving the crystal transver- 
sally in the beam tunes its birefringence and there- 
fore the OPO, which, therefore and unexpectedly, be- 
comes extremely sensitive to table vibrations trans- 
verse to its axis. We have confirmed this by observ- 
ing that the jitter's spectral signature matches that 
of the table vibration measured with accelerometers, 
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both sharing in particular sharp resonances at 70 and 
100 Hz. This produces noise on i'±, and Eq. © shows 
that the cavity length servo is ill-equipped to do much 
correction on because of the low tuning coefficient. 
As already mentioned, the temperature servo is too 
slow to be useful here, hence further narrowing of 
the beat note spectrum will require an electro-optic 
servo on the crystal so as to make use of the favorable 
tuning coefficient of Eq. {Tj) . 

Note, however, that the current level of perfor- 
mance is already outstanding and satisfactory for 
EPR and HLI experiments, the beat note being tun- 
able to exactly zero hertz (within the residual jitter 
range) by a simple temperature adjustment. 

3.3 Intensity control 

In order to determine how close to the threshold the 
OPO can stably be operated, which is crucial in order 
to determine whether EPR correlations can be ob- 
served, we have investigated the conversion efficiency 
given by [ST] 

^ ^ ^ Ps + P. ^ -(VlV-1), (9) 
^ Pp Pp ^ ' 

where N — Pp/Pth is the number of times above 
threshold (and should be as close to 1 as possible) 
and K is a, constant whose value is 2 for a standing 
wave cavity with a single-pass pump and up to 4 if 
the pump is reflected, the value of 4 corresponding 
to a ring resonator or to a double-pass pump with 
the right phase relationship with the signal and idler 
fields 6(j) — (j)p — 4>s — 4'i = — 7r/2 j^Hl- The value of K 
therefore depends on the differential phase shift im- 
parted to S(f) by reflection on the OPO output mirror 
and on air dispersion. We measured the OPO out- 
put power for different values of the pump power for 
the frequency-degenerate mode in stably locked con- 
ditions. Figure|Sldisplays a plot of p versus Pp, as well 
as a least-squares fit using Eq. I^. Since we do not 
know the mirror's phase shift, we leave K as a free 
parameter in the fit. The lowest point on the plot cor- 
responds to TV = 1.04, p = 5%, and Pout = 700 /xW 
per beam. We have in fact observed stable oscillation 
for N as low as 1.01, i.e. output powers of a few hun- 
dred microwatts per beam. This should indeed be 
in the regime where EPR correlations are expected 
[S]. The OPO output noise on each individual beam 
is mainly that of the pump laser, i.e. fairly low (but 
with relaxation oscillations) and is shot-noise limited 
above 1.5 MHz. As already seen on fig. 01 intensity- 
difference measurements are shot-noise limited above 
200 kHz. This speaks to the low level of classical in- 
tensity noise. Figure2|(a) also gives an idea of longer- 
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Figure 5: Measurement and fit of the conversion ef- 
ficiency of the OPO locked on the degenerate mode. 
p = ^{Vn — 1), where N = Pp/Pth- Fit parameters 
are Pth = 25.6(2) mW, K = 3.26(6), = 1-5 10^^^ 



term intensity fluctuations of the unlocked OPO in a 
single-detector measurement. 

The theoretical threshold is 12 mW ^T]. We mea- 
sured about 15 mW for the strongest cluster mode 
and 25.6(2) mW for the degenerate mode, which is 
not the strongest (center) cluster mode. This is due 
to a slight misalignment of the phase-matching angle 
of the crystal. Correcting it is quite doable but is 
always a delicate operation as it couples to the cavity 
alignment. We observe that the power level of the 
degenerate mode remains remarkably constant over 
time, even with the experiment (including the pump 
laser) turned off each night, testifying to the stability 
of the optical alignments. 

It is possible to achieve even lower powers with the 
same threshold, but the capture range of the servo 
decreases as the gain linewidth and this degrades the 
locking performance. We are currently implementing 
lower-noise detection techniques aiming at increasing 
the signal-to-noise ratio by one order of magnitude 
and at expanding the capture range of the servo. 
One should also note that, since the error signal is 
mixed with a considerable amount of vacuum fluctu- 
ations (99.99%), we only use a classical signal to lock. 
As a consequence, our servo cannot correct quantum 
fluctuations, be they phase difference noise or critical 
threshold fluctuations, which will allow us to observe 
them without interference from the locking electron- 
ics. 
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4 Conclusion 

We have experimentally demonstrated the stable op- 
eration of a frequency degenerate type-II OPO above 
threshold, with negligible polarization crosstalk, us- 
ing a noncritically phase-matched Sodium-doped 
KTP nonlinear crystal. The degenerate mode is ob- 
tained reproducibly without the need to realign the 
OPO cavity after the initial alignment, and its inten- 
sity is controlled as low as just a few percent above 
threshold. The twin beams exhibit 4 dB of number- 
difference squeezing at 200 kHz, which illustrates the 
low classical noise of our setup. Pushing this limit 
to even lower frequencies, below 1 kHz, is important 
from the perspective of applying HLI to the detection 
of gravitational waves. All stability and squeezing 
performances can be significantly improved on this 
starting point, which is already suitable for the gen- 
eration of bright EPR beams and HLI. 
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